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Abstract 

Background: Plasmodium vivax is one of the widespread human malarial parasites accounting for 75% of malaria 
epidemics. However, there is no baseline information about the status and nature of genetic variation of Plasmodium 
species circulating in various parts of Pakistan. The present study was aimed at observing the molecular epidemiology 
and genetic variation of Plasmodium vivax by analysing its merozoite surface protein-3a (msp-3a) and merozoite surface 
protein-3(3 (msp-3(l) genes, by using suballele, species-specific, combined nested PCR/RFLP detection techniques. 

Methods: A total of 230 blood samples from suspected subjects tested slide positive for vivax malaria were collected 
from Punjab, Sindh, Khyber Pakhtunkhwa, and Balochistan during the period May 2012 to December 2013. Combined 
nested PCR/RFLP technique was conducted using PvmspSa and Pvmsp-3[i genetic markers to detect extent of genetic 
variation in clinical isolates of P. vivax in the studied areas of Pakistan. 

Results: By PCR, P. vivax, 202/230 (87.82%), was found to be widely distributed in the studied areas. PCR/RFLP analysis 
showed a high range of allelic variations for both msp-3a and msp-3fi genetic markers of P. vivax, i.e., 21 alleles for 
msp-3a and 19 for msp-3(3. Statistically a significant difference (p < 0.05) was observed in the genetic diversity of the 
suballelic variants of msp-3a and msp-3/3 genes of P. vivax. 

Conclusion: It is concluded that P. vivax populations are highly polymorphic and diverse allelic variants of Pvmsp-3a 
and Pvmsp-3f3 are present in Pakistan. 

Keywords: Plasmodium vivax, PCR/RFLP, Genetic polymorphism, MSP3a and MSP3/3 genes 



Background 

Malaria is the major threat to public health and economic 
development in many nations [1]. It kills more than a mil- 
lion people a year, and approximately 40% of the world's 
populations live in malarious countries [2]. Plasmodium 
vivax is the most widespread species of human malaria 
parasites in the world and is endemic in many countries 
of Asia, Central and South America, the Middle East, and 
parts of Africa, with an estimated burden of 70-80 million 
cases annually [3,4]. In Pakistan P. vivax is the common 
malarial species, contributing to 70% of the malaria burden 
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[5]. Although P. vivax is highly prevalent, it has received 
scant scientific attention in Pakistan, thus, a paucity of 
baseline data exists on various aspects of P. vivax, such as 
population structure and drug resistance patterns. 

Genetic diversity is defined as the total number of 
heritable characteristics in the genetic make-up of the 
species [6]. Studies have shown that Plasmodium species 
show diversity in several parameters, including morphology, 
biochemistry, relapse patterns, symptoms, course and 
duration of infection, immunological responses, drug 
resistance, and transmissibility by anopheline vectors. 
However, the extent of genetic diversity observed depends 
on the transmission rates, immune pressure and natural 
selection of the parasite in an area [7] . Usually in endemic 
areas, where transmission rates are high, extensively di- 
verse variants are in circulation, posing a threat to public 
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health, since such variants have possibly acquired in- 
creased virulence and resistance to drugs in order to fit 
and survive [8] . 

Studies of the population structure of malaria parasites 
are important for understanding the evolution of para- 
site virulence and the role of parasite diversity in malaria 
transmission, and for designing control tools, including 
vaccines, as well as evaluating the impact of malaria con- 
trol measures [9,10]. 

PvmspSa and PvmspSfi have been used as markers in 
population genetic studies worldwide [11-14]. The best 
approach for detecting genetic diversity is to analyse 
more than one marker gene, because by doing so the 
probability of different clones sharing the same genotype 
accidentally, is substantially reduced [15]. Accordingly, it 
is suggested that analysis of both PvmspSa and Pvmsp3ji 
genes of P. vivax enables greater ability in identifying 
parasite haplotypes and detecting mixed strain infections 
[11,12]. Limited data regarding the genetic diversity of 
PvmspSa and PvmspSfi is available from one endemic 
area of Khyber Pakhtunkhwa province of Pakistan [16], it 
is therefore, hypothesized that new allelic variants may be 
present in other provinces of the country. This study will 
seek to elaborate on the genetic polymorphism of P. vivax 
malaria isolates and compare them with those reported 
from other parts of the world. 

Methods 

Sample collection 

This study was conducted with the approval of the Medical 
Superintendents of source hospitals and the ethics com- 
mittee of Kohat University of Science and Technology 
(KUST) Kohat Pakistan. A total of 230 blood samples 
were collected, after the informed consent from symp- 
tomatic patients, with the assessment of malaria control 
laboratories in the district headquarters hospitals of the 
four provinces: Punjab, Sindh, Khyber Pakhtunkhwa, 
and Balochistan during May 2012 to December 2013. 
Out of these 230 samples, 58 were collected from Punjab, 
54 from Sindh, 68 from Khyber Pakhtunkhwa, and 50 
from Balochistan. All the samples were screened for the 
detection of malarial parasites by microscope and PCR. 

DNA purification 

DNA of blood samples were extracted with the help of 
DNA extraction kit (GF1, Vivantus, USA). 

PCR amplification for Plasmodium vivax identification 

The nested PCR amplification procedures used by Snounou 
et al. and Singh et al. [17,18] were applied for amplification 
and detection of the 18S rRNA type A gene size 1.2 Kb, 
using species-specific primers rVTVl 5'-CGCTTCTAGC 
TTAATCCACATAACTGATAC-3 ' and rVIV2 5'-ACTTC 



CAAGCCGAAGCAAAGAAAGTCCTTA-3 ' for the de- 
tection of P. vivax. 

PCR of Pvmsp-3a genes 

By using PCR technique, the allelic variants of P. vivax 
merozoite surface protein-3oc were detected by standard 
protocol of Bruce et al. [11]. Primary PCR of PvmspSa 
was performed in 20 ml using 3 ml of DNA, using 
primers PI (5'-CAGCAGACACCATTTAAGG-3 ) and 
P2 (5 -CCGTTTGTTGATTAGTTGC-3'), while nested 
reactions were done in 20 ul with primers Nl (5'-GACC 
AGTGTGATACCATTAACC-3 ') and N2 (5 -ATACTG 
GTTCTTCGTCTTCA GG-3') using 3 ul of the primary 
PCR product. Taq DNA polymerase (2.6 units) was used 
and the following cyclic conditions were performed: 95°C 
for 3 min, 94°C for 30 sec, 58°C for 30 sec, 68°C for 
2.5 min (2, 30 times) and 72°C for 5 min. Two per cent 
agarose gel was stained with ethidium bromide in order to 
visualize the PCR separated product, under UV illumin- 
ation. Gel documentation system (Clearver Scientific, 
USA) was used for the photography. 

RFLP analysis of PvmspSa genes PCR products 

For the RFLP analysis of PvmspSa genes, Alul restriction 
enzyme (Fermentas, USA) was used. The total volume of 
the reaction medium was kept 20 ul having the following 
reagents Alul buffer (1.2 ul), PCR water (9.8 ul), Alul 
enzyme (1.0 ul) and product of mspSa gene (8.0 ul). 

PCR of PvmspSp genes 

The allelic variants of P. vivax merozoite's surface 
protein-3 were detected by PCR analysis, using standard 
protocol of Yang et al. [12]. For PvmspSfi, primary PCR 
was performed in 20 ul using 3 ul of DNA using primers 
PI (5 '-GTATTCTTCGCAACACTC-3 ') and P2 (5'-CT 
TCTGATGTTATTTCCAG-3 '). Nested reactions were 
done in 20 ml with primers Nl (5'-CGAGGGGCGAA 
ATTGTAAACC-3 ') and N2 (5 ' -GCTGCTTCTTTTGC 
AAAGG-3 ) using 1 ul of the primary PCR product as 
the template. Taq DNA polymerase (2.6 units) (Vivantis, 
USA) was used and the following cyclic conditions were 
performed: 94°C for 20 sec, 54°C for 30 sec, 68°C for 
2.5 min, 35 cycles. Agarose was used to visualize the 
PCR product under UV illumination. Gel documentation 
system (Clearver Scientific, USA) was used for the 
photography. 

RFLP analysis of PvmspSji PCR products 

The RFLP analysis of PvmspSfi gene was carried out, using 
Pest 1 restriction enzyme (Fermentas, USA) described by 
Yang et al. [12]. The total volume of the reaction medium 
was kept as 20 ul. The master mix for RFLP analysis was 
composed of Pest 1 Buffer 1.2 ul, PCR water 9.8 ul, Pest 1 
enzyme 1.0 ul and product of mspSji gene 8.0 ul. 



Khan ef al. Malaria Journal 2014, 13:355 
http://www.malariajournal.eom/content/13/1/355 



Page 3 of 7 



87.03 % 



88 % 




Punjab 
Balochistan 
Sindh 
KPK 



85.29 ,% 



| 91.37] % 

Figure 1 Overall prevalence of vivax malaria in Pakistan 
by PCR. 



Statistical analysis 

The data were statistically analysed by using SPSS 
(Version 20) with application of Kruskal Wallis test. 
The p value less than 0.05 was assumed to be statistically 
significant. 

Results 

Out of the total 230 samples, 202 samples were found 
positive for P. vivax both by microscopy and PCR. The 
distribution of the positive samples in the four prov- 
inces of Pakistan remained: 58/68 (85.29%) for Khyber 
Pakhtunkhwa, 53/58 (91.37%) for Punjab, 47/54 (87.03%) 
for Sindh and 44/50 (88%) for Balochistan (Figure 1). 

Analysis of pvmspSa gene by PCR /RFLP 

A sum of 130/202 samples of blood were successfully 
amplified for pvmsp-3a gene and were identified for 
different genotypes of the aforementioned gene. Three 



major allelic variants of different sizes, i.e., 2.2 kb (Type-A), 
2.0 kb (Type-B), 1.1 kb (Type-C) were detected. Besides 
these, an unusual band size of about 0.75 kb, named as 
Type-D and mixed genotypes designated as M were also 
observed (Figure 2). Out of the 130 amplified samples, the 
frequency of the different types of alleles remained 64/130 
(49.29%) for Type-A, 27/130 (20.76%) for Type-B, 27/130 
(20.76%) for Type-C, 9/130 (6.92 %) for Type-D, and 3/130 
(2.30 %) for mixed. Overall 12 suballelic variants (A1-A12) 
for Type A, 5 (B1-B5) for Type B and 4 (C1-C4) for 
Type-C were observed. The alleles A3, A9, Bl, B2, and 
C2 were frequently observed in all provinces of Pakistan 
(Table 1). 

Analysis of pvmsp-Z$ gene by PCR/RFLP 

In the current study 72/202 samples of blood had been 
successfully amplified for pvmsp-?>ji gene with a high 
frequency of polymorphism among their different allelic 
variants. The amplified DNA fragments for these alleles 
were classified according to their size variation. Based 
on variation in size of the fragments, the amplified sam- 
ples could be differentiated into three allele sizes: for Type 
A (2.0-2.5 kb), Type B (1.6-2.0 kb) and Type C (1.5 kb) 
(Figure 3). 

In all these amplified samples, a couple of samples 
reflect more than one band of different sizes. These mul- 
tiple bands show mixed infection, which revealed 2.77% 
of the allelic variants. Out of the total 72 amplified 
samples, 16 (22.22%) isolates for Type-A, 32 (44.44%) 
for Type-B and 22 (30.55%) for Type-C were observed. 
Among these isolates Type-B was found to have highest 
per cent frequency and Type-A was found the lowest in 
per cent frequency (Table 2). 

In this study pvmsp-?>(i gene was analysed using Pstl 
restriction enzyme and 7 suballelic variants, (A1-A7) for 



250bp 
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Figure 2 PCR/RFLP analysis of P. vivax msp-3a gene. V shows P. vivax (the amplified product is 120 bp). A, B, C and D indicate the main allele 
variants, while A, B, C and D with numbers indicate subtypes of the main allele variants. M is a 50 base pairs DNA marker and m is a 1 kilo base 
pairs DNA marker. 
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Table 1 Size, RFLP pattern and allelic frequency of Pvmsp-3a gene 

Size variation RFLP by A/u1 Alleles Data analysis n (%) frequency distribution of sub-allelic variants 









Punjab 
(n = 32) 


Sindh 
(n = 27) 


Khyber Pakhtunkhwa 
(n = 41) 


Balochistan 
(n = 30) 


Total 
(n = 1 30) 


A 


A1— A12 


1 


2 (6.25) 


2 (7.40) 


2 (4.87) 


2 (6.66) 


8 (6.15) 






2 


2 (6.25) 


1 (3.70) 


2 (4.87) 


1 (3.33) 


6 (4.61) 






3 


3 (9.37) 


2 (7.40) 


4 (9.75) 


2 (6.66) 


11 (8.46) 






4 


2 (6.25) 


1 (3.70) 


2 (4.87) 


1 (3.33) 


6 (4.61) 






5 


1 (3.12) 


1 (3.70) 


1 (2.43) 


1 (3.33) 


4 (3.07) 






6 


2 (6.25) 


1 (3.70) 


1 (2.43) 


1 (3.33) 


4 (3.07) 






7 


0 (0.00) 


1 (3.70) 


1 (2.43) 


1 (3.33) 


3 (2.30) 






8 


1 (3.12) 


0 (0.00) 


1 (2.43) 


1 (3.33) 


3 (2.30) 






9 


3 (9.37) 


3 (11.11) 


4 (9.75) 


3 (10) 


13 (10) 






10 


0 (0.00) 


0 (0.00) 


1 (2.43) 


1 (3.33) 


2 (1.53) 






11 


2 (6.25) 


0 (0.00) 


1 (2.43) 


1 (3.33) 


3 (2.30) 






12 


0 (0.00) 


0 (0.00) 


1 (2.43) 


0 (0.00) 


1 (0.76) 


B 


B1— BS 


13 


3 (9.37) 


3 (11.11) 


3 (7.31) 


2 (6.66) 


11 (8.46) 






14 


2 (6.25) 


3 (11.11) 


3 (7.31) 


3 (10) 


1 1 (8.46) 






15 


1 (3.12) 


1 (3.70) 


1 (2.43) 


1 (3.33) 


4 (3.07) 






16 


0 (0.00) 


0 (0.00) 


1 (2.43) 


0 (0.00) 


1 (0.76) 






1 7 


0 fO 00) 


0 (0 00) 


0 (0 00) 


0 (0 00) 


0 fO 00) 


C 


C1— C4 


18 


2 (6.25) 


1 (3.70) 


2 (4.87) 


1 (3.33) 


6 (4.61) 






19 


4(12.5) 


3 (11.11) 


4 (9.75) 


4 (13.33) 


15 (11.53) 






20 


2 (6.25) 


1 (3.70) 


2 (4.87) 


1 (3.33) 


6 (4.61) 






21 


0 (0.00) 


0 (0.00) 


0 (0.00) 


0 (0.00) 


0 (0.00) 


D 


D1 


22 


2 (6.25) 


2 (7.40) 


3 (7.31) 


2 (6.66) 


9 (6.92) 


M 


Ml 


23 


0 (0.00) 


1 (3.70) 


1 (2.43) 


1 (3.33) 


3 (2.30) 


Total frequency (%) 






32 (24.61) 


27 (20.76) 


41 (31.53) 


30 (23.07) 


130 (64.35) 


Mean ± SE 






1 .48 ± 0.80 


1 .30 ± 0.96 


1.1 7 ±0.92 


1.78 ±0.79 





Kruskal wallis test p < 0.05 

n = The number of subtypes of alleles that bear RFLP pattern. 
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Table 2 Size, RFLP pattern and allelic frequency of Pvmsp-3fi gene 



Size variation RFLP by Psf 7 


Alleles 




Data analysis n (%) frequency distribution of sub-allelic variants 








Punjab 


Sindh 


Khyber Pakhtunkhwa 


Balochistan 


Total 






(n = 21) 


(n = 20) 


(n = 17) 


(n = 14) 


(n = 72) 


A A1— A7 


A1 


1 (4.76) 


1 (5) 


0 (0.00) 


1 (7.14) 


3 (4.16) 




A2 


1 (4.76) 


1 (5) 


0 (0.00) 


0 (0.00) 


2 (2.77) 




A3 


1 (4.76) 


1 (5) 


1 (5.88) 


0 (0.00) 


3 (4.16) 




A4 


1 (4.76) 


2 (10) 


1 (5.88) 


1 (7.14) 


5 (6.94) 




A5 


1 (4.76) 


0 (0.00) 


1 (5.88) 


0 (0.00) 


2 (2.77) 




A6 


0 (0.00) 


0 (0.00) 


1 (5.88) 


0 (0.00) 


1 (1.38) 




A7 


- 








- 


B B1— B10 


Bl 


1 (4.76) 


1 (5) 


1 (5.88) 


1 (7.14) 


4 (5.55) 




B2 


2 (9.52) 


2 (10) 


2 (11.76) 


2 (14.28) 


8 (11.11) 




B3 


2 (9.52) 


2 (10) 


0 (0.00) 


1 (7.14) 


5 (6.94) 




B--1 


1 (4.76) 


1 (5) 


1 (5.88) 


0 (0.00) 


3 (4.16) 




B5 


1 (4.76) 


1 (5) 


1 (5.88) 


1 (7.14) 


4 (5.55) 




B6 


2 (9.52) 


3 (15) 


2 (11.76) 


0 (0.00) 


7 (9.72) 




B7 
B8 


0 (0.00) 


0 (0.00) 


0 (0.00) 


1 (7.14) 


1 (1.38) 




B9 


- 
_ 


- 
- 


- 
- 


- 
- 


- 
_ 




BIO 


- 


- 


- 


- 


- 


C CI— C4 


CI 


2 (9.52) 


2 (10) 


2 (11.76) 


1 (7.14) 


7 (9.72) 




C2 


2 (9.52) 


1 (5) 


1 (5.88) 


1 (7.14) 


5 (6.94) 




C3 


1 (4.76) 


1 (5) 


0 (0.00) 


1 (7.14) 


3 (4.16) 




C4 


2 (9.52) 


1 (5) 


1 (5.88) 


1 (7.14) 


5 (6.94) 




C5 


0 (0.00) 


0 (0.00) 


1 (5.88) 


1 (7.14) 


2 (2.77) 




C6 












M M1 


M1 


0 (0.00) 


0 (0.00) 


1 (5.88) 


1 (7.14) 


2 (2.77) 


Total frequency (%) 




21 (29.16) 


20 (27.77) 


17 (23.61) 


14 (19.44) 


72 (35.64) 


Mean ± SE 




1.11 ±0.17 


1.05 ±0.1 9 


0.89 ±0.1 5 


0.74 ±0.1 3 




Kruskal wallis test 








p < 0.05 






n = The number of subtypes of alleles that bear Psfl RFLP pattern. 



Type-A, 10 (B1-B10) for Type-B and 6 (C1-C6) for 
Type-C were observed. The suballeles B7-B10 were not 
found in the present study (Table 2). The different 
samples resolved through pstl restriction enzyme, show a 
high polymorphism in size and distribution among sam- 
ples (Figure 3). 

Discussion 

Pakistan is endemic for malaria and there are reports 
showing the distribution of various allelic variants of P. 
vivax. So, this study was aimed with the main objectives 
of identifying the new variants and to observe the epi- 
demiology of level of pvmspSa and pvmspSfi polymorph- 
ism of P. vivax in Pakistan. In the current study four 
distinct sizes of PCR products for PvmspSa were de- 
tected. Among these, three different major allelic variants 



of different sizes, i.e., 2.2 kb (Type-A), 2.0 kb (Type-B), 

1.1 kb (Type-C) along with unusual band size of almost 
0.75 kb or 300 bp name as Type-D and mixed genotypes 
designated as M, were observed in the study areas, 
while in Thailand [13], Afghanistan [19] and the Federally 
Administered Tribal Areas (FATA) of Pakistan [20], 
only three allele sizes (A, B and C) were obtained for 
this gene. 

The frequency of the different types of alleles 
remained 64/130 (49.29%) for Type-A, 27/130 (20.76%) 
for Type-B, 27/130 (20.76%) for Type-C, 9/130 (6.92%) 
for Type-D, and 3/130 (2.30%) for mixed, comparing 
with the study of Khatoon et al. where 82% (41/50) iso- 
lates having allele sizes 1.9 kb (Type-A), 6% (3/50) with 
the 1.5 kb fragment (Type-B), 8% (4/50) with the 

1.2 kb fragment (Type-C), 2% (1/50) with the 0.3 kb 
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fragment (Type-D) and one isolate (2%) exhibited mixed- 
strain infection [16]. 

In this study restriction digestion of the PvmspSa 
PCR product with Alul yielded 12 suballelic variants 
types: (A1-A12) for Type-A, 5 (B1-B5) for Type-B and 4 
(C1-C4) for Type-C were observed. The alleles A3 11 
(8.46%) and A9 13 (10%) for Type-A, allele Bl and B2 

11 (8.46%) for Type-B and C2 15 (11.53%) for Type-C 
were frequently observed in all provinces of Pakistan, 
while in the similar study carried out by Khatoon et al., 

12 different alleles designated as A1-A7, Bl, B2, CI, C2, 
and D, allele A3 being the most abundant 12/50 (24%) [16]. 

The PCR-RFLP patterns for P. vivax Pvmsp-3fi gene 
produced three categories Type-A (2.0-2.5 kb), Type-B 
(1.6-2.0 kb) and Type-C (1.5 kb). Similarly, in northwest 
Thailand, three distinct allele types (A, B and C (-0.65 kb)) 
were found [12]. In contrast, studies carried out in dis- 
trict Bannu by Khatoon et al. found two types of allele 
A (1.7-2.2 kb) and B (1.4-1.5 kb) of size polymorphisms 
for P. vivax Pvmsp-3ji gene [16]. Compared with P. vivax 
isolates from previous work conducted in Asia, the results 
of the current study further confirmed the existence of 
small geographic differentiation among P. vivax popula- 
tions. An earlier study in western Thailand showed that 
the B type of pvmsp-3fi is more abundant (60.4%) than 
other types, whereas in Chinese Bengbu and Guangxi 
samples, both A and B types were similarly prevalent [12]. 
According to Zhong et al., the Type-A allele was the most 
abundant in all four parasite populations (>57% more 
abundant than other types), i.e., Anhui, Hainan, Yunnan, 
and Myanmar [21]. In the current study, Type-B was 
found to have the highest per cent frequency and Type-A 
was found to have lowest per cent frequency. 

Data from PCR-RFLP of the Pvmsp-3a and Pvmsp-3f> 
loci showed that 5.02% of the Pakistani P. vivax isolates 
exhibited mixed-strain infections, which is comparable 
to sites in Bengbu (China) (5.6%) and is lower than 
observed in Thailand (20.5%) [12] and FATA of Pakistan 
(30%), which directly shares a border with Afghanistan 
[20]. In northern Iran [20] and Hongshuihe (China), no 
mixed genotypes were detected [12]. Apart from differences 
in transmission intensities in different areas, the observed 
variations in parasite types could be attributed to factors 
such as sampling bias, host immune selective pressure on 
particular types and/or spatio temporal changes in the 
availability of different mosquito species that can transmit 
specific parasite types in a particular area over different 
times or seasons, with areas that share the same mosquito 
types tending to have similar parasite types [20]. 

Conclusion 

This study indicates that P. vivax populations in four 
provinces of Pakistan: Punjab, Sindh, Khyber Pakhtunkhwa, 
and Balochistan, are highly diverse. Such heterogeneity 



might cause differences in parasite virulence, transmissi- 
bility and responses to chemotherapy, with important im- 
plications for malaria control measures in this populous 
region. 
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